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Abstract 

We summarize our recent results on the description of the evolution of Majorana neu- 
trinos in external fields. First we discuss the purely quantum method which involves the 
secondly quantized Weyl spinors. In frames of this formalism we find exact solutions of 
the classical Hamilton equations for Weyl fields in background matter, which then are used 
to canonically quantize the system. The analog of the effective Hamiltonian for the de- 
scription of spin-flavor oscillations of Majorana neutrinos in matter and a magnetic field 
is also obtained. Finally we discuss another approach for the treatment of Majorana neu- 
trinos in external fields, which is based on the relativistic quantum mechanics. The latter 
method does not require the quantization of the neutrino wave function. Nevertheless we 
demonstrate that the relativistic quantum mechanics approach also appropriately describes 
spin-flavor oscillations of relativistic Majorana neutrinos. 

1 Introduction 

It is known that the neutrino interaction with external fields can significantly influence the evo- 
lution of astrophysical [lj and cosmological media [2j . For example, the interaction of supernova 
neutrinos with a strong magnetic field of a protoneutron star (PNS) may cause a macroscopic 
effect like the kick of a pulsar [3]. The neutrino scattering off the matter of a rotating PNS 
may result in a spin-down of a star [4j. However the recent estimates [5] showed that the 
role of neutrinos in PNS spin-down is significantly less important than that in the possible 
explanation of great linear velocities of neutron stars. The neutrino interaction with primordial 
plasma is essential for the observed spectrum of the cosmic microwave background [6] and for 
the generation of large scale magnetic fields [7] . 

It was experimentally confirmed that neutrinos are massive particles and there is a mixing 
between different neutrino generations [8]. These neutrino properties result in the appearance 
of neutrino flavor oscillations [9] . It is also known that the presence of background matter can 
cause the resonance enhancement of neutrino oscillations [10] . Moreover the combined action 
of a background matter and an external magnetic field leads to the resonant transition like, 
v~ O ut (see, e.g., Ref. [H]), where the indexes ± denote different helicity states. Hence 
active neutrinos of the flavor a can be converted into sterile neutrinos of another flavor /?. 

Neutrinos can interact with an external electromagnetic field due to the presence of anoma- 
lous magnetic moments [32]. Note that the structure of the magnetic moments is completely 
different for Dirac and Majorana neutrinos (see, e.g., Ref. [13]). Despite the fact that nowadays 
there is no universally recognized confirmation of the nature of neutrinos [14] . in the present 
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work we shall suppose that neutrinos are Majorana particles. Note that in various scenarios 
for the generation of elementary particles masses it is predicted that neutrinos should acquire 
Majorana masses [15]. It should be also mentioned that numerous experimental attempts are 
made to investigate whether neutrinos are Dirac or Majorana particles |16j . 

In the present work we summarize our recent results [17U18] on the theory of Majorana neu- 
trinos oscillations in matter and electromagnetic fields. First, in Sec. [21 following the discussion 
of Ref. [20], we review the general properties of the neutrino mixing in vacuum. Then, in Sec. El 
we recall the basic equation for a massive Majorana neutrino propagating in background matter 
and electromagnetic field. In Sec. HI using the results of our recent works [THl[T9] we propose 
the classical field theory treatment of massive Weyl fields propagating in a background matter 
and interacting with an external electromagnetic field. Exact solutions of the wave equation 
for Weyl fields in a background matter are found in Sec. [5l In Sec. [6] we use these solutions 
to canonically quantize the Weyl fields. In Sec. in frames of our method we re-derive the 
effective Hamiltonian for the description of spin- flavor oscillations of Majorana neutrinos in 
matter and a magnetic field. Finally, in Sec. we consider the alternative approach for the 
description of spin-flavor oscillations of Majorana neutrinos based on the relativistic quantum 
mechanics and compare it with the results of Sec. In Sec. we discuss our results. 



2 Neutrino mixing in vacuum 

Experimentally it was found [21J that the number of active neutrinos is equal to three with the 
high level of accuracy. Nevertheless in various theoretical and phenomenological scenarios (see, 
e.g., Ref. [22]) some additional sterile neutrinos are discussed. Thus the general dynamics of the 
system of flavor neutrinos can be formulated in terms of the 3 + N s spinor fields, where N s > 
is the number of sterile neutrinos. It was experimentally confirmed [23] that active neutrinos 
correspond to the left-handed chiral projections, uh, A = e, /i,r, whereas sterile neutrinos 
should be expressed in terms of the right-handed chiral projections, z^, A = 1,. . . ,N S . Here 
^l,r = (1 =F7 5 )/2 x v. 

The system of wave equations for the spinors u^' will involve both Dirac and Majorana 
types of the mass terms. The simultaneous diagonalization of the mass terms can be performed 
with help of the matrix transformation [24] . 

N\ = Y, U X«Xa, (1) 

a=l 

where = (vj 1 , i/J, f^ 1 , [^fl c > ■•• > [^jy ] c ) is the multiplet of the neutrino fields, (U\ a ) is the 
unitary vacuum mixing matrix of the size (3 + N s ) x (3 + N s ), and Xa are the neutrino mass 
eigenstates. The charge conjugation is denned in the standard manner: v c = v^v* . 

The fields Xa are neither Dirac nor Majorana. However we may always construct Majorana 
fields of the the chiral projections of Xa as 

^a=M + ^(x«) c ]> (2) 

where x and yt are the phase factors having unit absolute values. Now we can see that the 
fields if) a satisfy the Majorana condition in an extended sense, 

tpa = x c if>a, (3) 
where x c = xx. In the following we shall suppose that x c = 1. 
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3 Electrodynamics of Majorana neutrinos in a background mat- 
ter 



The evolution of Majorana fields ip a , defined in Eq. ([2D, in a background matter under the 
influence of an external electromagnetic field is governed by the following wave equation: 

(iVfy - m a )Va - + <&7m7V& = 0, (4) 

where m a are the masses of the particles and = (i/2)[7^; 7^]_. Note that we will formulate 
the dynamics of the system Q in the mass eigenstates basis rather than in the flavor basis, as 
it is usually done when neutrino oscillations are considered, since, as we demonstrated in Sec. [21 
only in the mass eigenstates basis one can distinguish between Dirac and Majorana masses. 

The matrix of the neutrino interaction with a background matter, (g^ b ), introduced in Eq. Q 
in the mass eigenstates basis, is related to the analogous matrix in the flavor basis as follows: 

€b = E u lxfxx>Uy b , (5) 

A,A'=e,/i,r 

where the mixing matrix (U\ a ) is defined in Eq. (pTJ) . In general case the matrix (g^ b ) is hermitian. 
However we shall discuss the situation when the CP invariance is conserved. In this case the 
matrix (g^ b ) is symmetric. Despite a current attempt to detect CP violating terms in the 
neutrino sector [25], no definite results have been obtained yet. 

Note that in Eq. ([5]) we do not exclude the possibility of the existence of nonstandard neu- 
trino interactions, which correspond to the nondiagonal elements of the matrix (/vw) (see, e.g., 
Ref. [17J). In the situation when only the standard model interactions with matter are present, 
the matrix (f^y) is diagonal: f^ x , = 5\yf x , where f x are the effective potentials of interac- 
tions with background fermions. The zero components of these potentials, /P, are proportional 
to the effective matter density of non-moving and unpolarized background fermions, whereas 
the vector components, fA, are the linear combinations of the averaged matter velocity and 
the polarization. The explicit form of the effective potentials and the details of the statistical 
averaging can be found in Ref. [26] . 

Note that the vector term in the neutrino matter interaction ~ Oahlti^b 1S omitted in Eq. Q 
since it is washed out for Majorana neutrinos. The contribution of the axial-vector interaction 
with matter ~ g^^fi^ipb ^° ^ ne wave equation ([3D is twice the analogous contribution for Dirac 
particles since both neutrinos and antineutrinos equally interact with a background matter (see, 
e.g., Ref. (27J). 

Neutrinos can interact with the external electromagnetic field F^ v = (E, B) due to the 
presence of the anomalous magnetic moments (/i f,). It is know (see, e.g., Ref. [28]) that the 
matrix (/x a &) should be hermitian and pure imaginary, i.e. fi ab = —fJ-ba and fi* ab = —^ a b- Unlike 
the interaction with matter, which is generically defined for flavor neutrinos, the computation 
of the neutrino magnetic moments has to be done in the mass eigenstates basis. The detailed 
example of such a calculation is given in Ref. [29] , We shall discuss the situation when no 
admixture of sterile neutrinos is in the mass eigenstates ip a , i.e. when the eigenstates Z/? are 
absent. In this case the neutrino electric dipole moments are equal to zero |28j . 

In Ref. |30] is was rigorously proved that four component Majorana spinors are equivalent 
to two component Weyl spinors. Thus we may re-express the wave function ip a in Eq. ([5]) in 
two ways, 

<&> - (7") , » r # - (_^, ) . (6, 

Both representations in Eq. © satisfy the Majorana condition ([3]). Using Eq. ([6D we can rewrite 
Eq. © as 

fj a - {aV)r] a + m a a 2 rf a - Hab^(B - iE)a 2 r] b * + i{g Q ab + o-g ab )rj b = 0, (7) 
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or 

L + (<TV)£ a - m a a 2 C + Vab(T{B + iB)a 2 C b ~ Kdab ~ f*gab)6 = 0. (8) 

In the following we shall postulate these equations. Note that the analog of Eq. (J7J) was previ- 
ously derived in Ref. [31J. 



4 Classical field theory 

In Ref. |19j we demonstrated that the classical dynamics of a massive Weyl field in vacuum 
should be described only in frames of the Hamilton formalism. Generalizing the results of 
Ref. |19| to include the interactions with a background matter and an electromagnetic field we 
arrive to the following Hamiltonian (see also Ref. [18J): 

H = / d3r [E K ( ffV )% - (vlfi^X + m« [{rf a fa^ a + (K) T °2Va] } 

J a 

+ J2 {A*o6 ftMB - iE)(T 2 r/ 6 * + r]Ja 2 a{B + iE)<] 

a6 

- 1 [^I(9ab + 0-ga&)??6 " (^a) T (fafe + O-gaftX] } , (9) 

where 7r a are the canonical momenta conjugate to the "coordinates" r\ a . Using the aforemen- 
tioned properties of the matrices [}iab] and (5^) we find that the functional (J9j) is real as it 
should be for a classical Hamiltonian. 

Applying the field theory version of the canonical equations to the Hamiltonian H, 

5H 

Va =t — = (crV)r] a - m a a 2 if a + fi ab cr(B - \B)a 2 rf b - i(g® b + ag ab )r] b , (10) 
5H 

TTo = - 7— = (cr*V)iT a + m a a 2 ix* a - /i ab a 2 a(B + iE)^ + i(#° b + o-*g ab )ir b , (11) 

01Ja 

one can see that in Eq. (|10p we reproduce Eq. ([7|) for Weyl particles, which correspond to 
left-handed neutrinos, interacting with matter and an electromagnetic field. If we introduce the 
new variable ^ a = ia 2 TT a , we can show that Eq. ([IT]) is equivalent to Eq. ((HJ) for right-handed 
neutrinos. 



5 Exact solution of the wave equation 

In the following we shall suppose that the background matter in average is at rest and unpo- 
larized, i.e. g ab = 0. This approximation is valid in almost all realistic cases (see the detailed 
discussion in Ref. |18j). Let us decompose the Hamiltonian ([9]) into two terms H = Hq + H{ nt . 
The former term, Hq, contains the vacuum Hamiltonian as well as the matter term diagonal in 
neutrino types, 

H = I d 3 r £ K T [(<rV) - iflgjifc - (v* a ) T [(<rV) ~ i<7°J< 

^ a 

+ m a [{rfcfwa + «) T a 2 77 a ] } . (12) 
The latter term in this decomposition, 

= [ d 3 r £ {Hob [vrJ<x(B - iEy 2% * + vIwiB + iE)vr 6 *] -ig° ab [vrjr/, - (^*) T <] }, (13) 

J a^b 

has the nondiagonal matter interaction and the interaction with an electromagnetic field which 
is nondiagonal by definition. 
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Analogously to Eqs. (jlOp and (jlip we define the reduced Hamilton equations which contain 
only the Hamiltonian Hq: r}a = SHq/Stt^ and tt^ = — 5Hq/5i]^\ Using the results of 
Refs. |18y31j we can find the solutions of these equations in the form, 



a„ W-e 



d 3 p 

(2vr) 3 / 2 

(at)*W-e iEit + 



m n 



E+ + 



w+e 



-iEZt 



D ipr 



(a-y W+ e iE ^ 



E a + |p| +g\ 



d 0) (r,0 



d 



(2tt) 3 / 2 1 

(b-y w+ e iE ^ 



b^w+e 



-iE+t 



+ 



Et + |p| - 9 c a 



Ea + \p\+g° a 

-(b+Tw-e 



-b„ W-e 



-ipr 



= ipr 



iEtt 



-ipr 



(14) 



where we introduce the new variable £, a °^ = ^2^a°\ w± are the helicity amplitudes defined in 
Ref. [32], and 

(15) 



4° = V^ + (IpI-C0 2 - 



is the energy of a Weyl field [27J, £ = ±1 is the particle helicity. To derive Eqs. (|14p and (|15|) 
we suppose that the external field g® a is spatially constant. 



6 Quantization 

Now we can carry out the canonical quantization of the Weyl fields r] a and ir a . For this purpose 
we suggest that the expansion coefficients a^(p) and h„ (p) in Eq. (fT4"|) are operators. Note 
that the operators in the expansion of ij a and £ a are independent since in Sec. H] we showed that 
these fields evolve independently. 
Using the following relation: 

a±(p)(£± + |p| T 9° aa ) = £(p)(|p| T 9 ° aa ), (16) 

and suggesting that the operators a„ (p) obey the anti-commutation properties, 

{a±(k);[a±( P )]*}+ = ^ 3 (k-p), (17) 

with all the rest of the anticommutators being equal to zero, we can cast Eq. (|12p into the form, 

H = d 3 p£[£~ (a~)*a~ + E% (a^)* 0+] + divergent terms, (18) 

J a 

which shows that the total energy of a massive Weyl field is a sum of energies corresponding to 
elementary oscillators of positive and negative helicities. Note that the details of the derivation 
of Eq. CED can be found in Ref. [E] 

Using the results Ref. |19j we can also quantize the total momentum of a Weyl field defined 

as 



(19) 

With help of Eqs. (|14p . (I16p . and (|17p we rewrite Eq. (|19p in the following form: 

P = / d 3 p£ p[(a-)*a- + (a+)*a+] + divergent terms, (20) 

which has the similar structure as Eq. (|18j) . 
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It is interesting to mention that a massive Weyl field in vacuum can be quantized in the 
two independent ways (see Ref. [19]) because of the degeneracy of the neutrino energy levels: 
E~ = E+ = \J + |p| 2 . On the contrary in matter only one of the possibilities for the 
quantization gives the correct result for the total energy (|18p since the energy levels are no 
longer degenerate, cf. Eq. (fl~5|) . 



7 Nondiagonal interaction with matter and a magnetic field 

To quantize the Hamiltonian H[ Qt we shall use the forward scattering approximation. It means 
that one has to account for only the terms conserving the number of particles [33J. Using 
Eqs. (fT4"|) . (fl~6j) . and (fTTj) we rewrite Eq. (fT3]) in the form, 

H int = [ d 3 p^[M afe (a-)*a fe -e i ^ + M+(a+)*a+e i ^ 

+ F^a-ya+e^ + F+ia+y^e^}, (21) 

where «5± = E± - E± and o± = E± - Eg. 

The general form of the coefficients Mi and F~t is quite cumbersome. Therefore we describe 
the dynamics of the system in the ultrarelativistic approximation, |k| 3> max(m a , g® a ), where k 
is the initial momentum of neutrinos. Moreover we discuss the simplest case of the two neutrino 
eigenstates, a = 1, 2, and suppose that E = 0, since it is difficult to create a large scale electric 
field. In this limit we get for the coefficients M h ?» ^fg a b and F^ — // a b|B| sin-^kB) where 
??kB is the angle between the vectors k and B. The explicit form of these coefficients for the 
arbitrary neutrino momentum can be found in Ref. |18j . 

Now we define the neutrino density matrix as 

S 3 (p - k)p AB (k) = (a* B (p)a A (k)), (22) 

where A = (£, a) is a composite index and (...) is the statistical averaging over the neutrino 
ensemble. Applying the quantum Liouville equation for the description of the density matrix 
evolution, 

ip=[p,H int ], (23) 

we can rewrite it as ip = [H,p], where 



is the effective quantum mechanical Hamiltonian. 

To study the evolution of our system we use Eq. (|14p . where the wave functions already 
contain time dependent exponential factors. The energies in Eq. (|14p correspond to the total 
diagonal Hamiltonian (|12p . cf. Eqs. (|15p and (|18p . which contains both the mass term and the 
diagonal interaction with a background matter rather than only a kinetic term as in Ref. |34j . 
Thus our treatment is similar to the Dirac picture of the quantum theory. That is why in 
Eq. (|23p it is sufficient to commute the density matrix only with H int rather than with the total 
Hamiltonian H = Hq + H{ nt . 

To eliminate the time dependence in the effective Hamiltonian (|24|) we make the transfor- 
mation of the density matrix |31j . 

Pqm = UpU\ U = diagje-^+^i)*, e^*-*^*, eT^-^i)*, e^+^t}, (25) 

where $ = 5m 2 /4|k| is the phase of vacuum oscillations and 5m 2 = m^ — m 2 , is the mass squared 
difference. 
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The evolution of the transformed density matrix can be represented as ip qm = [H qm , Pqm], 
where the new effective Hamiltonian has the form, 



9°2i 



/ $ + g°n 9°i2 

-^ + 522 -^2i|B|sin^ kB 
-/zi 2 |B|sin^ kB $-911 
\-/x 2 i|B| sini9 kB 



-9°2i 



-yui 2 |B| sin i? kB \ 


-g°i2 

-$-9°22 J 



(26) 



Recalling the properties of the magnetic moments matrix for Majorana neutrinos, [112 = i/U and 
H21 = — !/■*> with fi being a real number, we can see that Eq. (|26p reproduces the well known 
quantum mechanical Hamiltonian for spin- flavor oscillations of Majorana neutrinos in matter 
and a magnetic field [TT] . 



8 Relativistic quantum mechanics 

In Sec. [7] we showed that the dynamics of neutrino spin-flavor oscillations in matter and a 
magnetic field can be described in frames of the purely quantum approach based on the quantum 
Liouville equation (|23|) for the neutrino density matrix (|22|) . There is, however, an alternative 
formalism for the treatment of neutrino oscillations which involves the relativistic quantum 
mechanics method |17y31| . This approach is based on the exact solutions of wave equations for 
massive neutrinos in an external field, which were obtained in Sec. [5j 

The main idea of the relativistic quantum mechanics method is the following: for the given 
initial wave functions of flavor neutrinos one should find the wave functions at subsequent 
moments of time. Note that this approach is similar to the relativistic wave packets description 
of neutrinos oscillations [35] . However we study the propagation of neutrino wave packets which 
exactly account for external fields. If one defines some known initial distribution of active 
neutrinos, then we can calculate the transition and survival probabilities for other neutrino 
flavors. Note that using Eq. ([1]) we can always express the initial condition for flavor neutrinos 
via that for the neutrino mass eigenstates, i.e. we shall suppose that r) a (r,t = 0) = r{i (r) are 
the given functions. 

In our analysis we shall assume that the expansion coefficients in the decomposition of the 
wave functions (|14p are c-numbers rather than operators as it was supposed in Sec. (6) Thus, 
since we do not quantize the Weyl fields, we no longer need the additional degree of freedom 
represented by the spinor £. Moreover we may change the normalization of the spinors at our 
convenience. 

We will look for a general solution of Eq. (J7|) in the following form (see Eq. (|14p and Ref. [31): 

+ f[a+(p,t)]W^ + F ™° o KM]*w + ^*) e-'A, (27) 
V E a + |p| +<, / J 

where the expansion coefficients a„ (p, t) are the functions of time to account for the nondiagonal 
terms in Eq. ([7]). The energy levels and the helicity amplitudes w± are defined in Sec. 

We choose the plane wave initial wave functions of massive neutrinos, rj^ (r) = A a e lkr , where 
A a is the amplitude of the initial neutrino field distribution. If we are not interested in the next- 
to- leading effects in neutrino oscillations, we may discuss the ultrarelativistic approximation 
from the very beginning: |k| 3> max(m a , g^a)- As in Sec. [7] we suggest that E = as well as 
adopt the following geometry: B = Be z and k = |k| x (sin^kBj 0, cos $ k B)- 
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Inserting the ansatz (|27p in Eq. ([7]) and after the straightforward calculations we get the 
system of the ordinary differential equations for a„ (p, i). It is convenient to rewrite it in a con- 
ventional form as the Schrodinger-like equation, i^f = H^, introducing the quantum mechanical 
"wave function" of neutrinos, Vl/ T = (a^ , , • • • , ajj" , , ■ ■ • ), and the effective Hamiltonian H, 
which has the same form as in Eq. (|24p . As in Sec. [7J here we shall discuss the situation of 
two neutrino eigenstates, a = 1,2. Then we make the transformation, ^ qm = U*&, where the 
matrix U is defined in Eq. (|25|) . The evolution of the transformed quantum mechanical "wave 
function" $ qm is governed by the effective Hamiltonian H qm which coincides with that given 
in Eq. (|26|) . Therefore we obtain that the relativistic quantum mechanics method is consistent 
with the standard quantum mechanical description of neutrino spin-flavor oscillations in matter 
and a magnetic field. 

9 Conclusion 

In the present work we have summarized our recent results on the theory of spin-flavor oscilla- 
tions of Major ana neutrinos in matter and a magnetic field. After the brief introduction to the 
theory of neutrino mixing in vacuum in Sec. [2] and the discussion of the Majorana neutrinos 
electrodynamics in Sec. [3j we have mainly considered two approaches: the purely quantum 
method, Sees. 0]|71 and the relativistic quantum mechanics approach, Sec. El 

The former method is based on the canonical quantization of Weyl fields in presence of a 
background matter. The procedure of the canonical quantization requires the construction of 
the classical field theory Hamiltonian [36J , in which all the dynamical observables should then be 
replaced with operators. Previously it was claimed [37] that massive Weyl fields are essentially 
quantum objects to be expressed via anticommuting operators. Nevertheless, in Sec.Ul we have 
constructed the classical field theory of massive Weyl fields in a background matter and an 
electromagnetic field, which is based on the Hamilton formalism. Note that one does not doubt 
that generically our world is quantum. However there are numerous processes which may be 
also described within the classical physics (see Ref. [55] for many interesting examples). 

In Sec. [5]we have found the exact solution of the wave equation for a Weyl field in presence of 
a background matter. This solution has been used in Sec. [6] to canonically quantize the system. 
Then, in Sec. [7J we have completed the quantization of the system by taking into account the 
nondiagonal interaction with matter and a magnetic field using the density matrix formalism. 
For ultrarelativistic particles we have re-derived the effective Hamiltonian for the description 
of spin-flavor oscillations of Majorana neutrinos in matter and a magnetic field. 

Finally, in Sec. [HJ we have considered the alternative approach for the description of the 
evolution of Majorana neutrinos which is based on the relativistic quantum mechanics. In 
frames of this method we do not quantize the neutrino wave functions. Instead of operators 
in the expansion of the fields, we deal with c-number functions. These functions evolve in 
time in such a way to satisfy the given initial fields distribution. We have demonstrated that 
for ultrarelativistic neutrinos this method also gives the appropriate description of spin-flavor 
oscillations of Majorana neutrinos in external fields. 
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